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ABSTRACT 

Recent ultraviolet spectra from the Hubble Space Telescope (HST) and Far 
Ultraviolet Spectroscopic Explorer (FUSE) satellites demonstrate that UV line 
and continuum fluxes observed from Mira B are increasing back towards the levels 
that the International Ultraviolet Explorer observed in 1979-1980 and 1990-1995, 
after having been found in a very low state by HST and FUSE in 1999-2001. The 
UV emission is associated with accretion of material onto Mira B from Mira A's 
massive wind, so the variability is presumably due to variations in accretion rate. 
From wind absorption features, we estimate a Mira B mass loss rate of 2.5 x 10~ 12 
M Q yr _1 , indicating that Mira B's wind has increased in strength along with the 
accretion rate. The UV variability may be associated with a previously reported 
14-year periodicity in Mira B's optical emission. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual 
(o Ceti) — stars: winds, outflows — ultraviolet: stars 



1. INTRODUCTION 

Mira AB is the nearest symbiotic system at a distance of only ~ 130 pc, (Perryman et 
al. 1997), and the only symbiotic binary in which the components of the system have been 
resolved at wavelengths ranging from X-rays to radio (Karovska et al. 1997, 2005; Matthews 
& Karovska 2006). The primary (Mira, o Ceti, HD 14386) is a very large pulsating AGB 
star (photospheric diameter of ~ 500 R©), with a massive, cool wind of ~ 10~ 7 M© yr _1 
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(Bowers & Knapp 1988). The companion (Mira B, VZ Ceti) is accreting material from 
Mira A, but the exact nature of the accretor is uncertain (Jura & Helfand 1984); it has 
usually been assumed to be a white dwarf (Joy 1954; Karovska et al. 2005). The current 
angular separation between the components is ~ 0.6" (Karovska et al. 1997). The orbit 
of the binary is poorly known, but Prieur et al. (2002) estimate a long orbital period of 
~ 500 years, and this orbit suggests that the distance from Mira A to Mira B is close to the 
apparent separation in the plane of the sky, about 70 AU. This interacting binary, being one 
of the few accretion systems in which the components of the system are resolvable, provides 
a unique opportunity to study the spectral energy distributions of the components of an 
interacting binary individually, as well as being a laboratory for testing accretion models. 

In 1995, the Faint Object Camera (FOC) instrument on the Hubble Space Telescope 
(HST) resolved the Mira system for the first time spatially and spectrally at UV and optical 
wavelengths (Karovska et al. 1997). Since then, the system has been monitored from UV 
to radio wavelengths. In the past few years, we have been witnessing dramatic changes in 
Mira B, especially at UV wavelengths (Karovska et al. 1997; Wood et al. 2001). Spectra from 
the Space Telescope Imaging Spectrograph (STIS) instrument on HST imply a substantial 
decrease in accretion luminosity in the late 1990s, compared with fluxes observed at earlier 
times by the International Ultraviolet Explorer (IUE) and HST/FOC. Furthermore, a forest 
of Lya-fluoresced H 2 emission lines dominated the HST spectra of Mira B in 1999, despite not 
being seen at all in any previous IUE or HST spectra (Wood et al. 2001, 2002). Observations 
from the Far Ultraviolet Spectroscopic Explorer (FUSE) demonstrated that the low accretion 
luminosity and prominent H 2 emission persisted into 2001 (Wood & Karovska 2004). On 
2003 December 6, Chandra observations detected an unprecedented soft X-ray outburst from 
Mira A (Karovska et al. 2005), further evidence of remarkable dynamicism in the Mira AB 
system. The outburst was also detected in contemporaneous Ha observations and is likely 
associated with a magnetic flare followed by a mass ejection or jet-like activity (Karovska et 
al. 2005, Karovska 2006, in prep.). 

We here present recent UV observations of Mira B from HST and FUSE, obtained fol- 
lowing the X-ray outburst. We compare these data with past observations by these satellites 
and by IUE. We discuss the long term variations of Mira B and a possible correlation with 
long term optical variability. 

2. UV SPECTROSCOPY OF MIRA B FROM 1979-2001 

Numerous UV spectra of Mira AB were obtained in 1979-1980 and 1990-1995 by IUE. 
This data set includes 32 spectra taken in IUE's short-wavelength, low- resolution mode (i.e., 
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SW-LO spectra, covering ~ 1150 - 1950 A). These far-UV (FUV) spectra are dominated 
by broad, hot lines such as C II A1335, Si IV A1400, C IV A1550, Si III] A1892, and C III] 
A1908, which are produced by material accreting onto Mira B from Mira A's massive wind. 
The few high-resolution spectra of this wavelength region (i.e., SW-HI spectra) are very 
noisy, but they suggest very broad widths for these lines, consistent with emission from a 
rapidly rotating accretion disk (Reimers & Cassatella 1985). Unlike IUE, HST can resolve 
the binary, and UV spectra from HST confirm that the broad lines are from Mira B, and 
that Mira A contributes very little to the IUE spectra below 2000 A (Karovska et al. 1997; 
Wood et al. 2001). 

Only a modest level of UV variability is present within the large IUE SW-LO data 
set, with line and continuum fluxes varying by about a factor of 2. However, when HST 
observed Mira B in 1999 August using STIS, the character of Mira B's FUV spectrum had 
changed dramatically. The UV continuum and all of the aforementioned UV emission lines 
had weakened by at least a factor of 20, and the FUV spectrum was now dominated by many 
narrow H2 lines fluoresced by H I Lja (Wood et al. 2001), which are not detectable at all in 
the IUE spectra. The accretion onto Mira B results in a warm, fast wind from the star that 
is most clearly apparent in broad absorption features in the Mg II h & k lines at 2800 A. 
Given the presence of accretion, this wind is envisioned to be of a bipolar nature. Analysis 
of the Mg II absorption demonstrates that the substantial drop in UV flux from Mira B seen 
by HST is accompanied by a corresponding drop in its mass loss rate (Wood et al. 2002). 

The low UV fluxes observed by HST in 1999 are most naturally explained by a sub- 
stantial decrease in the accretion rate onto Mira B, which also led to a weaker wind from 
the star. The sudden prominence of the H 2 lines in the HST spectra suggests that unlike 
other UV lines, the H I Lya line responsible for fluorescing these lines had not decreased 
in strength relative to the IUE era. Studying Lyct in IUE spectra is very difficult due to 
geocoronal contamination, but Wood et al. (2002) still uncovered evidence to support the 
contention that Lyct fluxes are similar in the HST and IUE spectra. 

A 2001 November observation from FUSE found Mira B to be still in a very low state, 
with C III A1175 fluxes comparable to those in the 1999 HST/STIS spectrum (Wood & 
Karovska 2004). Like the STIS data, this FUSE spectrum also shows numerous narrow H 2 
lines. Ground-based observations have suggested that Mira B's optical light curve may have 
a period of 14 years, with the last published minimum being in 1971 (Joy 1954; Yamashita 
& Maehara 1977). Wood & Karovska (2004) noted that the the numerous IUE SW-LO 
observations are all near predicted maxima of this cycle, while the 1999 HST and 2001 
FUSE observations are near a predicted minimum. Thus, it is possible that the observed 
long-term UV variability could be correlated with previously detected optical variations. If 
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the optical periodicity is real and if a connection between the optical and UV light curves 
truly exists, Mira B's UV fluxes should increase towards a 2006-2007 maximum predicted by 
the 14-year optical periodicity. The new 2004 HST and FUSE spectra presented here allow 
us to test this prediction, and to see if the character of Mira B's UV spectrum has changed 
once again. 

3. RECENT HST AND FUSE OBSERVATIONS 

3.1. Spectral Line Variability 

Mira B was observed by HST with the STIS instrument on 2004 February 16. A near- 
UV (NUV) spectrum from 2274-3119 A was taken using the moderate resolution E230M 
grating, with an exposure time of 600 s, and Mira B's FUV spectrum from 1140-1729 A 
was observed using the E140M grating, with an exposure time of 2100 s. Both spectra were 
taken through the 0.2" x 0.2" aperture. 

Mira was also observed by FUSE on 2004 August 31. This was an 8603 s exposure 
through the LWRS aperture, covering the 905-1187 A spectral range. This represents only 
a fraction of the requested 30 ksec exposure, so the new FUSE spectrum is much noisier 
than the 2001 FUSE observation (Wood k Karovska 2004). A C III A1175 flux can still be 
measured from the data and some of the strongest H 2 lines can still be discerned. 

Figures 1 and 2 compare UV spectral regions observed by HST in 2004 with the old 
1999 spectra. Five of the Figure 1 panels focus on spectral regions containing FUV emission 
lines of H I, O I, C II, Si IV, and C IV. Some of these lines are blended with H 2 emission, 
so the locations of strong, narrow H 2 lines are also indicated in the figure. The last panel 
shows the NUV spectrum in the region of the Fe II UV1 multiplet lines. The line locations 
indicated in Figure 1 are with respect to Mira's rest frame, where the rest wavelengths are 
taken from Abgrall et al. (1993) and Morton (2003), and we assume a radial velocity of 56 
km s" 1 for the Mira system (Bowers & Knapp 1988; Planesas et al. 1990; Josselin et al. 
2000). 

Figures 1 and 2 show that UV fluxes from Mira B have increased dramatically since 
1999, with the O I, C II, Si IV, C IV, and Mg II lines being much stronger than they were. 
The notable exceptions to this rule, the H I Lya line and the H 2 lines fluoresced by Lya, 
will be discussed further in §3.3. The UV continuum is also much stronger in 2004 (see the 
Fe II panel in Fig. 1). In Figure 3, we plot measured fluxes as a function of time for many 
different lines, including measurements from IUE, HST, and FUSE. For C III A1175, O I 
A1300, C II A1335, Si IV A1400, and C IV A1550, the IUE fluxes are measured from the 32 
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IUE SW-LO observations mentioned in §2. Examples of these spectra have been previously 
presented by Reimers & Cassatella (1985) and Wood et al. (2001). Note that the fluxes in 
Figure 3 are the total flux of these multiplets (e.g., for Si IV the flux includes both the A1393 
and A1402 lines), given that the individual components of the multiplets are not resolved in 
the IUE SW-LO spectra. 

Blended H 2 lines are removed from the data before measuring fluxes from the HST and 
FUSE data. This cannot be done for the low resolution IUE spectra, but thanks to the higher 
line fluxes seen by IUE, H 2 is not as strong a contaminant in the IUE data, so the correction 
is not as important (Wood et al. 2001). The C III A1175 line is the only line observable 
by all three satellites, though in the HST/STIS spectra and the 2004 FUSE spectrum the 
C III line is observed at very low S/N. In the higher quality 2001 FUSE spectrum, the H 2 
lines add 1.6 x 10~ 14 ergs cm -2 s _1 to the C III feature (Wood & Karovska 2004), so we 
simply subtract this flux from each C III flux measurement to correct for the strong H 2 
contamination in C III. 

The IUE's inability to resolve the Mira binary is generally not a problem since Mira A 
does not contribute significantly to most UV emission lines. However, Mira A can contribute 
to Mira AB's Mg II lines at certain pulsation phases (Wood & Karovska 2000). Fortunately, 
in long-wavelength, high-resolution (LW-HI) IUE spectra, Mira A's narrow Mg II lines are 
mostly contained within the saturated core of Mira B's Mg II wind absorption. Therefore, 
the Mira A Mg II line can be removed before measuring a Mg II flux for Mira B. The IUE 
Mg II fluxes plotted in Figure 3 are all from LW-HI spectra. We do not use LW-LO spectra 
because of potential Mira A contamination. The Mg II fluxes in Figure 3 include both the 
h and k lines. 

Figure 3 shows the dramatic decrease in UV line fluxes observed in 1999-2001. The 
highest temperature line, C IV A1550, exhibits a particularly large decrease of nearly a 
factor of 100. The fluxes increase significantly in 2004, but they are still lower than fluxes 
observed by IUE in 1979-1980 and 1990-1995. As discussed in §2, a 14-year periodicity has 
been detected in Mira B's optical light curve (Joy 1954; Yamashita & Maehara 1977), and 
dotted lines in Figure 3 indicate the phase and period of this light curve. The increase in 
UV flux that we see in our 2004 spectra is roughly consistent with the predictions of the 
optical light curve, assuming that the UV emission follows the same periodic behavior as that 
inferred for the optical light. It therefore remains a plausible hypothesis that Mira B's long- 
term UV variability is associated with the previously detected periodic optical variability. 
However, truly verifying this connection would require UV observations over at least another 
cycle. 
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3.2. Wind Absorption 

The blue sides of all the resonance lines in Figures 1 and 2 are suppressed by absorption 
from Mira B's wind, with the exception of C IV (see below). The wind absorption is more 
readily apparent in the new 2004 data than in the old 1999 spectra. This is in part due to 
much better S/N in the new data thanks to the higher fluxes, but it is also due to the wind 
absorption being intrinsically stronger, which is particularly evident for the Fe II and Mg II 
lines. The mass loss rate of Mira B's wind can be estimated most easily from the Mg II wind 
absorption profiles, which have the highest S/N. 

From analyses of the Mg II k line at 2796 A, we have previously estimated a mass loss 
rate of M = 5 x 10~ 13 M Q yr" 1 from the 1999 HST spectrum, and M = 1 x 10~ n M yr" 1 
from a typical IUE spectrum (Wood et al. 2002). Using the exact same methodology, we 
have fitted a wind absorption profile to the 2004 Mg II k line in Figure 2. The observed 
absorption profile is not smooth, suggesting a more complicated and structured wind than 
we can fit with such a simple model. Nevertheless, the mass loss rate inferred by the fit 
should still be a reasonable estimate. This fit suggests a mass loss rate of M = 2.5 x 10~ 12 
M Q yr _1 and a wind terminal velocity of \4o = 450 km s -1 . For comparison, the absorption 
predicted by the wind absorption profile from 1999 is also shown in Figure 2, assuming 
M = 5 x 10~ 13 M Q yr" 1 and = 250 km s^ 1 . The wind absorption is significantly stronger 
in 2004, meaning that Mira B's wind has clearly strengthened, but the mass loss rate is still 
not as high as the M — 1 x 10~ n M yr _1 value estimated for the IUE era. The mass loss 
rate variations are consistent with the accretion rate variations inferred from the UV fluxes 
(see Fig. 3). Clearly, Mira B's accretion rate is highly correlated with the strength of its 
wind. 

The only atomic lines in Figure 1 that do not definitively possess wind absorption are 
the C IV lines. Small decreases in flux blueward of the center of the two C IV lines may 
indicate the presence of some C IV wind absorption, but this absorption is not nearly as 
prominent as the saturated absorption seen for the other lines. The presence of strong wind 
absorption for the Si IV A 1392 line, and the relative lack of absorption for C IV implies 
a maximum temperature for the wind between T = 6.0 x 10 4 K (Si IV's line formation 
temperature) and T — 1.0 x 10 5 K (C IV's line formation temperature). 

3.3. The Anomalous Lya Line 

The one clearly anomalous line in Figure 1 is H I Lya. Although the 2004 Lyct flux 
is slightly higher than in 1999, this line does not show the dramatic flux increase that the 
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other emission lines show. The singular lack of variability for this line was also noted by 
Wood et al. (2002), who inferred that the Lya flux during the IUE era cannot have been 
much stronger than in 1999, unlike all the other lines. The lack of variability for Lya means 
that the H 2 lines that are fluoresced by Lya are also relatively constant, as shown by the H 2 
lines blended with C II and Si IV in Figure 1. 

It is difficult to understand why the Lya; line behaves so differently from the continuum 
and all other emission lines, including lines formed at similar temperatures (O I, C II, and 
Mg II). Does this mean that the Lya emission is coming from a different place than the 
other lines? If so, why? Another possibility is that the Lya flux actually does vary along 
with the other lines, but the extremely high opacity of Mira B's wind to Lya moderates the 
variation. Wind absorption explains why the observed Lya emission is entirely redshifted 
relative to line center (see Fig. 1). Analogous absorption is also seen for many of the other 
lines, but wind opacity will be orders of magnitude higher for Lya since hydrogen is much 
more abundant than other atomic species. 

Given that increases in UV flux are clearly accompanied by increases in wind strength, 
perhaps the enhanced wind obscuration largely cancels the flux increase for Lya. This 
hypothesis would require that Lya photons are actually destroyed rather than just scattered 
by the wind. Perhaps this destruction could be accomplished by scattering into Mira B's 
accretion disk, or extinction by dust entrained within Mira B's wind. 

4. CAUSES OF LONG-TERM VARIABILITY 

The fundamental cause of Mira B's substantial variability remains uncertain. Variable 
mass loss from Mira A could have caused the changes in the UV flux of Mira B observed in 
the past 10 years, signaling that the system is undergoing important transformations. These 
transformations may have set the stage for the remarkable X-ray outburst from Mira A that 
was observed by Chandra on 2003 December 6 (Karovska et al. 2005). 

If the variability is associated with the possible 14-year optical periodicity, the reason 
for the periodicity is a mystery. The UV lines in particular are believed to originate from 
accretion processes on Mira B, so a periodicity in these lines would suggest a periodicity for 
the accretion rate. Accretion instabilities could easily lead to variability, but it is not clear 
why they would be periodic. Perhaps the gravitational influence of an undetected stellar 
or planetary companion could lead to a corresponding periodicity in the accretion rate, 
perhaps through wind focusing effects. Hints of a third body in the Mira AB system have 
been previously reported in the literature (Baize 1980; Karovska 1992; Karovska et al. 1993). 
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Since Mira A's wind is the source of the accreting material, a periodicity in Mira A's mass 
loss rate could in principle also lead to periodicity in the accretion rate. No such periodicity 
has yet been found in the light curve of Mira compiled by AAVSO, but searches for very low 
amplitude, long-term periodicity are ongoing (M. Templeton, private communication). 

We hope to explore these possibilities further in the future, as we reach the expected 
maximum in 2006-2007 and over the next 14-year cycle, using different observational tech- 
niques including optical/IR interferometry. However, future UV high resolution spectra are 
very much in doubt due to the failure of the STIS instrument in 2004 August, and the reac- 
tion wheel and gyro failures on FUSE that currently limit observations to high declinations. 
Thus, new UV instrumentation may be necessary to continue spectroscopic observations in 
the ultraviolet. 

This work was supported by NASA grant NAS8-39073 to the Smithsonian Astrophys- 
ical Observatory and NASA grant NNG05GC14G to the University of Colorado. MK is a 
member of the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical 
Observatory under contract to NASA NAS8-39073. 
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Fig. I. — A comparison of UV spectra of Mira B taken in 1999 (dotted line) and 2004 (solid 
line). The spectra are smoothed slightly for the sake of appearance. Dashed lines indicate 
the rest wavelengths (in Mira's rest frame) of the line or lines indicated by the panel labels 
(e.g., H I Lyct, O I, C II, Si IV, C IV, and Fe II). Dot-dashed lines indicate the locations of 
strong, narrow H 2 lines. 
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Fig. 2. — A comparison of Mira B's Mg II k line profile from 1999 (dotted histogram) and 
2004 (solid histogram), plotted on a velocity scale in Mira B's rest frame. A wind absorption 
profile has been fitted to the 2004 data (thick solid line), which suggests a mass loss rate 
of M = 2.5 x 10~ 12 M Q yr" 1 and a terminal velocity of = 450 km s _1 . Also shown is 
the absorption predicted by a fit to the 1999 wind absorption (dot-dashed line), assuming 
M = 5 x 10~ 13 M yr' 1 and = 250 km s" 1 (Wood & Karovska 2004). This greatly 
underpredicts the amount of absorption observed in 2004, illustrating the increase in wind 
strength from 1999 to 2004. 
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Fig. 3. — Fluxes of various UV emission lines as a function of time, based on spectra from IUE 
(triangles), HST (boxes), and FUSE (diamonds). The fluxes are normalized to the average 
value observed by IUE. The lines and normalization factors (in units of ergs cm -2 s _1 ) are 
C III A1175: 4.1 x 1(T 13 , O I A1300: 8.1 x 1(T 13 , C II A1335: 6.5 x 1(T 13 , Si IV A1400: 
5.5 x 1(T 13 , C IV A1550: 1.7 x 1(T 12 , and Mg II A2800: 1.3 x 10" 11 . The dotted lines 
are schematic representations of the 14 year periodicity of Mira B suggested by optical 
observations (Joy 1954; Yamashita & Maehara 1977). 



